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Abstract 

In this note, we present a new generator for Deeply Virtual Compton Scattering processes. 
This generator is based on formalism of Generalized Partons Distributions evolved at Next Leading 
Order (NLO). In the following we give a brief description of this formalism and we explain the 
main features of the generator for the elastic reactions, as well as for proton dissociation. 



1 Introduction 




Figure 1: The contributing diagrams for the reaction ep e^p : (a) diagram for the Deeply Virtual 
Compton Scattering process ; (b) and (c) diagrams for Bethe-Heitler process. 

Deeply Virtual Compton Scattering (DVCS) is the exclusive production of a real photon in 
diffractive lepton — hadron interactions, l + N ^ l + j + N, as shown in Fig. ^(a)- This process is 
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calculable in perturbative QCD (pQCD), when the virtuality, = — g^, of the exchanged photon 
is large and it interferes with the purely electro-magnetic Bethe-Heitler (BH) reaction presenting 
the same final state (Fig. [U (b) and (c)). 

The first measurements of the DVCS process at high W"^ = {q + P)"^ D3 12] and its beam-spin 
asymetry in polarised ep scattering at low 0] have recently become available. A proposal 
for dedicated studies at the COMPASS experiment (in /ip collisions) is also under review. 

The DVCS reaction can be regarded as the elastic scattering of the virtual photon off the 
proton via a colourless exchange. The pQCD calculations assume that the exchange involves 
two partons, having different longitudinal and transverse momenta, in a colourless configuration. 
These unequal momenta are a consequence of the mass difference between the incoming virtual 
photon and the outgoing real photon. The DVCS cross section depends, therefore, on the Gen- 
eralised Parton Distributions (GPD), which carry information about parton correlations inside a 
nucleon. 

First predictions for DVCS cross sections were based on LO calculations [Sj. These calculations 
explicitely involve the variable R defined as the ratio of the imaginary part for the amplitudes for 
DIS and DVCS processes : 

Im{A{YP 7*P)) 
Im{A{'~f*p '-fp)) 

A value of i? ~ 0.55 has been shown to be in good agreement with HERA results. Previous 
Monte-Carlo (MC) for DVCS generation are based on this approach [HIIZI- We refer to this pre- 
diction as FFS (Frankfurt-Freund-Strikman) in the following [5]. 

However, even if this effective LO prediction reproduces correctly the experimental data, it 
is not sufficient as it does not provide a direct insight of the rich information present in GPDs. 
The MC described in this note has been developped to allow experimental measurements to be 
compared with GPD models and to study asymetries. 

GPDs have been studied extensively in recent years |HJ [HI 1101 • These distributions are not 
only the basic, non-perturbative ingredient in hard exclusive processes such as DVCS or exclusive 
vector meson production, but they are also generalizations of the well known Parton Distribution 
Functions (PDFs) from inclusive reactions. 
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In the approach developped in reference [H^, the scattering amphtude is simply given by the 
convolution of a hard scattering coefficient computable to all orders in perturbation theory with 
one type of GPD carrying the non-perturbative information. The formalism for the NLO QCD 
evolution equations for GPDs can be found in reference 

Higher twists contributions to the DVCS amplitude have been calculated in reference |Sj. They 
consist in terms of order 0{mj^/Q), 0{y/^/Q) and 0{\qcd/Q) with m^r the proton mass and 
t = {P — P'Y the momentum transfer to the outgoing proton. They have been shown to be sizeable 
at low within the kinematic domain relevant for HERMES, CLAS or COMPASS experiments. 

2 Structure of the MC program 

The core of the program consits of several fortran routines provided by A. Freund. The calculations 
of cross sections from GPDs are done in two steps : 

1. The GPDs H, E, E are evolved at NLO by an independent code jlSj which provides tables 
for the real and imaginary parts of the so called Compton Form Factor (CFFs) : at LO, 
they are just the convolution of GPDs by the coefficient functions 8J. For example 

m, Q\ t) = J2 [' — 7^— ± ^ -owx, Q^ t)dx 

with Cj the fractional quark charge and C, = xb/{'2 — xb) the skewing variable. There are 
different tables for each GPD and for LO/NLO approximations, as well as for the twists- 
3 corrections. The GPDs parameterisations at the initial scale (in Q^) are described in 
reference 0. 

In the DGLAP domain, they consist in distributions based on the forward CTEQ6 param- 
eterisations ; with no external skeewing added in a profile function. In each table, there 
are 48 bins in xb = extending from 10~^ till 0.7 and for each of them 40 bins in 
from 1 GeV^ till 10^ GeV^ (regularly spaced in InQ^). We obtain the values for the real and 
imaginary parts of the CFFs considered after a 2D spline of the xb — grid. 

2. From these CFFs, the cross sections for DVCS and BH-DVCS interference are calculated 
according to formulae of reference (see sectional). The BH cross section is also calculated 
according to standard expressions. During the iterative integration step, the code provides 
an output with the estimated cross-section and the accuracy of the integration. It is always 
important to check the proper convergence of the integrals and that the final accuracy is 
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small (below 1%), which means that no problem have occured in the integration over the 
kinematic range considered. 

The squared amplitudes are first integrated within the kinematic domain defined by the user. 
Events are then generated according to the differential cross sections. These two steps make use 
of the BASES/SPRING package For each accepted event, the kinematics of the final state 
particles is stored into a PAW ntuple. The relevant kinematic formuae are detailed in section 
The generated events are stored in the standard LUJETS common of the PYTHIA [I5j program. 

3 Cross section formulae and kinematics 

In the MC, we are calculating the five-fold cross section for the DVCS process 

da (y^XBV 
dxBdyd\t\d(f)d(p IGn'^Q'^y^l + 

in which the amplitude T is the sum of the DVCS and Bethe-Heitler (BH) amplitudes : T = 
'^vcs + '^H- This cross section depends on the Bjorken variable xb, the squared momentum 
transfer t = {P — P'Y, the lepton energy fraction y = P ■ q/ P ■ k, and, in general, two azimuthal 
angles. We use the convention e = 2xb^. 

In equation (jT)), = 0Ar — 0/ is the angle between the lepton and hadron scattering planes 
and (fi = $ — (pN is the difference of the azimuthal angle <P of the transverse part of the nucleon 
polarisation vector S, i.e., S± = (0, cos^, sin<P, 0), and the azimuthal angle (pN of the recoiled 
hadron. Our frame is rotated with respect to the laboratory one in such a way that the virtual 
photon four-momentum has no transverse components, see Fig. El For the kinematics we choose 
the following convention : the z-component of the virtual photon momentum is negative and x- 
component of the incoming lepton is positive with k = {E, E sin 9i, 0, E cos 9i), q = (go, 0, 0, — l^sl). 
Other vectors are P = (M, 0, 0, 0) and P' = {E2, IP2I cos^sin 6'Ar, IP2I sin0sin6'jv, IP2I cos6'7v)- 

3.1 Cross sections 

The BH amplitude is real (to the lowest order in the QED fine structure constant) and is pa- 
rameterised in terms of electromagnetic form factors, which we assume to be known from other 
measurements. 

According to reference 0, the BH, DVCS and interference terms of equation (^J, namely 
I^BHp, I^DVcsP and J = Tbvcs'^bh + '^vcs'^bh can be written : 



J 
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Figure 2: The kinematics of the leptoproduction in the target rest frame. In this figure taken from 
reference fc, k' , P2 and q2 are respectively the four-momenta of the initial lepton, the scattered lepton, 
the outgoing proton and the real photon. The z-direction is chosen counter-along the three-momentum 
of the incoming virtual photon. The lepton three- momenta form the lepton scattering plane, while the 
recoiled proton and outgoing real photon define the hadron scattering plane. In this reference system the 
azimuthal angle of the scattered lepton is i;^; = 0, while the azimuthal angle between the lepton plane and 
the recoiled proton momentum \s (pN = (j). When the hadron is transversely polarised (in this reference 
frame) S±_ = (0, cos sin(^, 0), the angle between the polarisation vector and the scattered hadron is 
denoted as if = ^ — (pN- 
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where the + (— ) sign in the interference stands for the negatively (positively) charged lepton 
beam. In this expression, the Fourier coefficients q, Si are fonction of CFFs (and then GPDs) and 
1/Vi and 1/7^2 are the lepton BH propagators. 
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We have mentioned above that the real and imaginary parts of the CFFs (convolution of the 
GPDs by a hard coefficient function) are given in tables of xb — Q'^, thus the t dependence is 
assumed to be factorised and has to be defined by the MC user. We have included an option in 
the steering which allows to consider a global exponential behaviour of the cross section 

da/dt oc exp{B{Q^) t) 

in which the t slope B{Q'^) can be chosen as a linear function of \nQ^. Alternatively, the t 
dependence can also be given by the Pauli-Dirac form factors for each GPD fiavor. 

3.2 Event kinematics 

The reconstruction of energies and angles of the scattered positron, photon and proton is carried 
out for each event from x, Q^, t and (p. 

The kinematics of the outgoing lepton (and then of the virtual photon) is trivially obtained in 
the laboratory frame. The kinematics of the outgoing nucleon 

P' = {E2, IP2I cos sin 6* AT, IP2I sin0sin6'7v, I-P2I cos6'Ar) 
is most easilly obtained in the the frame of Fig. ^ where the target nucleon is at rest. We use 

E -m ' IP I ^t{t-Am%) 
J^2 — rriN — \-t^2\ — 7: 

and 

"^''" = 2^5^^^ 

in which the virtual photon components can be expressed as 



_ t + Q^xb 

2mN 2mNXB\ Q 



After a Lorentz transformation, we finally get all variables (energies and angles of incoming 
and outgoing particles) in the laboratory frame. In the ntuple produced during the generation, 
these variables are given in both frames. 
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4 Integration and Generation 



Integration and generation are realized with the BASES/SPRING package ^1]. By BASES, 
probabihty distributions are calculated by integrating the differential cross-section over the phase 
space and saved in a file which is then used in the generation step. By SPRING, events are 
generated by the MC method according to the above distributions. 

5 Proton dissociation 

It has been shown that the proton dissociation is a non negligible contribution (10% to 20%) in 
the measured DVCS samples for HI and ZEUS experiments E] . Thus, it is important to include 
this process in the MC. 

In case of proton dissociation the DVCS process ep — >• 61^7 leads to a state of mass My and 
the generation of the My spectrum follows the parameterisation exposed in ^Hl) as it is done in 
the DIFFVM MC [T7]. The main hypothesis is that can be factorised with the elastic DVCS 
cross-section, with : 

da _ 
dM^ M2(i+^) 

In the continuum region {My > 3.6 GeV^), /(M-p-) = 1 leading to a global My dependence in 
l/M^*^^"*"^) for proton dissociation. In the resonance region (My < 3.6 GeV^), /(My) is the result 
of a fit for the proton diffractive dissociation on deuterium (at fixed t) ^H]- It is important to 
note that this treatment is only acceptable for DVCS at large W {W > 20 GeV). Thus, it must 
not be applied in the BH mode and in the low W domain, where the stucture of resonances is 
more complex. 

6 Examples 

6.1 Elastic DVCS process 

To illustrate the ouput of the DVCS MC, we first present some examples in the elastic case 
{ep — s> eprf). In Fig. |21we present the generated kinematic variables in the range 10~^ < xbj < 0.1, 
4 < < 100 GeV^ and t > —1 GeV^. The generation is done at NLO for different t dependences 
of the GPDs as explained in the caption of Fig. 01 For the same generation. Fig. |3] represents the 
spectra for energies and polar angles of the outgoing particles. The normalization is done to the 
number of events. 
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Figure 3: Predictions for the kinematic variables at NLO of the DVCS MC for different behaviours in t 
: the fuh histogram stands for a global exponential t dependence of the DVCS cross section in * with 
B = 7 GeV~^ ; the full line represents also an exponential t dependence with a slope of S = 4 GeV~^ 
and the dotted line has been obtained when considering the Dirac-Pauli form factors for each GPDs. 
Histograms are normalized to the number of events. 

In Fig. 13 we show the coplanarity (absolute difference of the azimuthal angles of the scattered 
lepton and real photon) of the elastic DVCS process (again with the three different t behaviours 
considered since Fig. We observe the broadening of this distribution when the t slope becomes 
larger. Indeed, for small momentum transfer the ej system is balanced and the coplanarity is 
close to 180", whereas for larger momentum transfer the ej system becomes unbalanced and a 
deviation form 180° is observed. And, of course, the fraction of events with a large momentum 
transfer is enhanced when B becomes smaller. 

As mentioned in section ^ previous measurements (and then acceptance corrections) at low xb 
[1 , 2\ are based on the FFS approximation for the DVCS cross-section [5j. On Fig. IHl we compare 
the FFS (for R=0.55) and NLO predictions of the DVCS MC in the same generated kinematic 
range : 10~^ < xbj < 0.1, 4 < < 100 GeV^ and t > — 1 GeV^. For this comparison we 
only consider a global exponential t dependence with B = 7 GeV~^. Histograms are normalized 
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Figure 4: With the same convention as in Fig. El (different t dependences), we represent the energies 
and angles of the scattered (outgoing) lepton Elq and Olq-, for the outgoing proton Epo and 6po and 
for the real photon Erg and Org- 

to the number of events and we notice a good agreement in shapes. However the cross section 
calculated in the NLO case is calculated to be about 15% higher than the LO FFS one (for the 
GPDs parameterisations included in the code). Regarding the good agreement in shapes, one can 
wonder whether it is useful to extend the MC from FFS LO to NLO GPDs calculations. However, 
it is important to note at this stage that GPDs at NLO become essential to predict asymetries as 
explained later in this note. Also, as mentioned previously, with GPDs we can consider different 
t dependences for the different flavors, which is obviously not possible in the FFS approximation. 
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Figure 5: Coplanarity : absolute difference of the azimuthal angles of the scattered lepton and real 
photon, for the three cases mentioned in Fig. |31 
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Figure 6: Comparisons of the predictions based on the LO calculations 5^ with a value of i? = 0.55 
(black points) and of the NLO calculations (full histogram). In both cases, we have considered a global 
exponential t dependence with B = 7 GeV~^. Histograms are normalized to the number of events. 
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6.2 DVCS vs BH cross-sections 



In Fig. [7|we present another illustration of the MC : we compare the DVCS with the BH process 
when the lepton is scattered backward {6io > 160°), the real photon is scattered with 6j.g < 160° 
and the azimuthal angle is integrated over. This configuration is the one analysed in HI and 
ZEUS experiments to extract the DVCS cross-sections Histograms (Fig. [7j) are normalized 

to luminosities and the interference between DVCS and BH is calculated to be negligible : as the 
azimuthal angle is integrated over, the constant term in the expression for the interference (see 
formula 0]) is the only (negligible) contribution. We notice the different behaviours in W, Eio, 
which allow a good separation power at low W (W < 120 GeV) and large Eio {Eio > 15 GeV) 
and thus allow the measurement of the DVCS cross-section. 




W (GeV) E,^ (GeV) 




Figure 7: Spectra of DVCS (full line) compared to BH (full histogram) processes for W, energies of the 
outgoing lepton and the real photon and the coplanarity. The configuration considered for these plots is 
the one with the lepton scattered backward {Oi^ > 160°) and the real photon scattered with Org < 160°. 
Histograms are normalized to luminosities as calculated by the MC. 
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In Fig. IHlwe compare DVCS and BH when both lepton and photon are scattered backward 
{(^io,rg > 160°). Histograms are normahzed to luminosities. We notice that in this case the BH 
process is dominating, the DVCS representing about 6% of the sample and the interference about 
0.5%. 




Figure 8: Spectra for DVCS (full line) compared to BH (full histogram) processes for W, energies of 
the outgoing lepton and the real photon and the coplanarity. The configuration considered for these 
plots is the one with the lepton and the real photon scattered backward {Oio^rg > 160°). Histograms are 
normalized to luminosities as calculated by the MC. In this configuration, it is interesting to notice that 
the interference (dotted line) is small but non negligible (about 0.5%) of the BH sample. 
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6.3 Proton dissociation DVCS process 

As mentioned in section |31 in case of proton dissociation the DVCS process ep — > eY'j leads to a 
state of mass My. We present on Fig. El the generated My for two values of the t slope (slope of 
the global exponential t dependence). As mentioned in sectional we notice the My behaviour in 

da /(Ml-) 
dM^ M2(i+^) 

with the resonance region at low My (My < 3.6 GeV^) and the shape at larger values. 




tGeV 
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Figure 9: Example of predictions of the DVCS MC for proton dissociation. With two values of B (slope 
of the global exponential t dependence) -.3 = 1 GeV^^ (black line) and B = 2 GeV"^ (full histogram), 
we represent the t spectrum and the generated mass of the proton dissociated system My (in log). 

Away from the resonance region, the system Y is treated as a quark-diquark (q-qq) system. 
Its hadronisation is performed by PYTHIA. We have added an option compared to the treatment 
done in the DIFFVM MC. In DIFFVM, it is assumed that the proton splits into a q-qq system, 
so that the quark couples to the pomeron leaving a diquark remnant. Due the possible spin states 
of this diquark system (in a singlet or triplet wave function), different probabilities are assigned 
to the different configurations of the q-qq system. However, for a DVCS process, the coupling of 
the valence quark can also be done to a quark : see for example the LO handbag diagram. Thus, 
the probabilities for the q-qq system describing the proton have to be modified by taking into the 
account the electric charge coupling. We have introduced a real parameter in the steering which 
allows to define the type of probabilities considered : a = if we only consider the probabilities 
defined by spin states Pq , a = I'm. case of spin states and electric coupling Pi, and < a < 1 for 
P„ = (1 - a)Po + aPi- 
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7 Steering cards 

EXP 1 : 1/0 to use de logarithmic/linear binning in the kinematic range for Q"^ and xb 
SEED 2345671 : generation seed value 

FIXED FALSE : FALSE/TRUE to select collider/fixed-target modes 
ELEP -27.55 : beam momentum (in GeV/c) of the lepton beam 

ETARG 920. : in collider mode, beam momentum (in GeV/c) of the proton beam or, in fixed 
target mode, mass of the target 

LCHAR +1. : lepton charge in units of e in case TINTIN is set to FALSE (the code is using the 
grids of real and imaginary parts for CFFs) 

LPOL 0. : polarisation of lepton, +1 in direction of movement, -1 against direction of movement, 
for unpolarised 

TPOL 0. : target polarisation, +1(-1) = polarised along (opposite) probe/target beam direction, 
for unpolarised 

ZTAR 1. : charge of the target 

ATAR 1 . : atomic number of the target 

SPIN 1. : spin of target nucleon type (in units of /i/2) 

IRAD : 1 for QED radiative corrections to be apphed, otherwise (only the ISR are included 
in the code) 

lELAS 1 : 1 to run the code for the elastic case only, to run in the proton dissociation mode 
only. Note that wc can not run both cases together and, as mentioned above, the proton 
dissociation treatment implemented in the code can only be used for a pure DVCS process 

IRFRA 1 : in case of proton dissociation for the resonances domain, IRFRA=0 if the user wants 
to use DIFFVM's routines for the decays of resonances, 1 if this user wants to decay via 
PYTHIA 

PROSPLIT 1. : treatment of dissociated proton in the continuum domain, PROSPLIT = a (see 
above) steers how the proton is splitted into quark-diquark 
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EPSM 0.08 : for proton dissocation, EPSM = e with da/dM^ oc 1/My^^+'^ 

TINTIN FALSE : TRUE to run in the LO EES approximation, EALSE to use the grids for 
CEEs 

BTIN 2. : slope of the exponential t dependence (in the EES approximation) 
RTIN 0.55 : value for the variable R = ^i^^Aly^p^-yp)) (^^^ section [TJ 

E2QCD TRUE : in the EES approximation (Jdvcs oc F2/R'^, this parameter determines the 
parameterisation considered for F2 : TRUE/EALSE respectively for the HI QCD fit or 
ALLM results 

DIPOLE EALSE : if the parameter TINTIN is TRUE, it's possible to run the code using the 
diople model formula (in this case DIPOLE must be set to TRUE) 

IGEN : generation mode with BASES/SPRING : for the grid calculation and the generation 
and 4 for the generation only (which requires the file bases. data) 

NGEN 10 : dummy 

NPRINT 10 : number of events for which the PHYTHIA output is printed 
NCALL 10000 : parameter for BASES, number of sampling points per iteration 
ITMXl 10 : parameter for BASES, number of iterations for the grid defining step 
ITMX2 10 : parameter for BASES, number of iterations for the integration step 
IDEBUG : debug flag 

NXGRID 58 : Number of X points in the CEEs grids 
NQGRID 40 : Number of points in the CEEs grids 

IPRO 2 : process to generate : 1,2,3,4 or 5 for BH, DVCS, interference BH-DVCS, charge 
asymetry or single spin asymetry respectively. 

lORD 2 : order for the calculations of the GPDs amplitudes : 1/2 for LO/NLO 

XMIN l.Oe-4 : lower bound of for the kinematic domain of the generation. If the MC user 
run the code with the use of CEEs grids, this value must be larger than the lower bound of 
the grids (10"'^ in our case) 
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XMAX l.Oe-1 : upper bound of xb 
QMIN 3.0 : lower bound for (jn GeV^) 
QMAX 300.0 : upper bound for (jn QeV^) 
TMIN 0.0 : upper bound for t ^ {p' - pf (in GeV^) 
TMAX -1.5 : lower bound for t ^ {p' - pf (in GeV^) 
MYMIN 1.13 : lower bound on My in proton dissociation case 
MYMAX 30. : upper bound on My 

NSET 2 : integration over PHI = 0:2 for integration, 1 otherwise 

PHI 0. : in case of NSET=1, azimuthal angle between lepton and final state scattering planes 
PPHI 0. : in case of NSET=1, angle between lepton plane and transverse polarisation vector 
THETA 0. : angle between longitudinal and transverse polarisation 
TWIST3 FALSE : Twists 3 calculations (TRUE) or not (FALSE) 

ITFORM : in case of the use of CFFs grids, if ITFORM=0, the t dependence is in exp(6t) with 
h = BQCST + BQSLOPE logiQ^QOSQ) ; if ITF0RM=2, the Pauli-Dirac form factors are 
used for the t dependence ; finally, if ITF0RM=4, the t dependence is in exp{bt) for Xb < 
XODEF and Pauli-Dirac form factors are considered otherwise. Note that in case of the LO 
FFS approximation the t dependence is set by the parameter BTIN (see above). 

BQCST 7.0 : see ITFORM 

BQSLOPE 0.0 : see ITFORM 

BGCST 7.0 : see ITFORM 

BGSLOPE 0.0 : see ITFORM 

QOSQ 2.0 : see ITFORM 

XODEF O.leO : see ITFORM 

THLMIN 00.0 : lower bound on Oi^pton (in deg.) 
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THLMAX 180.0 : upper bound on Oiepton (in deg.) 

ELMIN 10.0 : lower value of E^is (in GeV) 

THGMIN 00.0 : lower bound on (in deg.) 

THGMAX 180.0 : upper bound on 9^ (in deg.) 

EIMIN 0.00001 : minimum value for the ISR photon (in GeV) 

YMAX 0.4 : maximum value for y = Q^/^sx) 

M12MIN 1. : minimum value for the reconstructed mass of the outgoing lepton and the real 
photon 

ELMB 0. : minimum value of the outgoing lepton energy (in GeV) 
EGMB 0. : minimum value of the real photon energy (in GeV) 

8 Summary 

MILOU is a new generator for DVCS based on the formalism of Generalized Partons distributions. 
This MC has been dcvclopped to allow experimental measurements to be compared with GPD 
models (at LO or NLO) and to study asymctrics. The generation of BH processes and interference 
between DVCS and BH are also aviable. In addition, in case of pure DVCS we have included the 
possiblility to study proton dissociation. 
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